Bulk Fermi surface and momentum density in heavily doped La2 xSi"a;Cu04 using high 
resolution Compton scattering and positron annihilation spectroscopies 
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We have observed the bulk Fermi surface (FS) in an overdoped (a;=0.3) single crystal of 
La2-iSr2,Cu04 by using Compton scattering. A two-dimensional (2D) momentum density recon- 
struction from measured Compton profiles yields a clear FS signature in the third Brillouin zone 
along [100]. The quantitative agreement between density functional theory (DFT) calculations 
and momentum density experiment suggests that Fermi-liquid physics is restored in the overdoped 
regime. In particular the predicted FS topology is found to be in good accord with the correspond- 
ing experimental data. We find similar quantitative agreement between the measured 2D angular 
correlation of positron annihilation radiation (2D-ACAR) spectra and the DFT based computations. 
However, 2D-ACAR does not give such a clear signature of the FS in the extended momentum space 
in either the theory or the experiment. 

PACS numbers: 74.72-h, 78.70.Ck, 78.70.Bj, 71. 10. Ay 
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I. INTRODUCTION 

High temperature cuprate superconductorgl^ such as 
La2-a;Sra;Cu04 (LSCO) show interesting Mott insulator 
behavior— 1^ at low doping x while at high doping these 
materials display Fermi liquid propertiesJ^i^ This fas- 
cinating phase diagram is complicated by the presence 
of various inhomogeneities such as stripe ordering^ and 
possible clustering of dopants.® Many models of the over- 
doped metallic phase have been put forward, some gen- 
eralizing the concept of the Fermi liquid, while others at- 
tempt to use Hubbard or t-J models. In this context, it is 
important to investigate the electron momentum density 
of these materials, which, when compared to predictions 
of theoretical models, can give an indication of their cor- 
rectness or applicability. 

Insight into the evolution of electronic states can be 
obtained via Fermi surface (FS) measurements in various 
doping regimes. Experimental FS work on the cuprates 
has to date been limited mainly to angle-resolved pho- 
toemission (ARPES),^'^ quantum oscillations (QO)jS 
and scanning tunneling spectroscopies (STS)»^ However, 
ARPES^^ and STS^^ are surface sensitive probes. Al- 
though QO measurements probe the bulk, they provide 
only FS areas without giving information on the loca- 
tion of the FSs in momentum space. Furthermore QO 
require long mean free paths and large magnetic fields 
which could alter the ground state. For underdoped sam- 



ples there is evidence of FSs distorted from Local Density 
Approximation (LDA) based band structure predictions, 
involving arcs or pockets, but as doping is increased it 
appears that a large, LDA-like FS becomes manifest as 
the pseudogap collapses near optimal doping. For an 
overdoped Tl-cuprate, the FS has been found to be large 
and closed around (tt, 7r)fi- suggesting that the Van Hove 
singularity (vHS) still lies below the Fermi energy Ep in 
good agreement with LDA calculations!^ 

These considerations provide motivation for deploying 
genuinely bulk sensitive spectroscopies for FS measure- 
ments in the cuprates. Two such spectroscopies, which 
have been used extensively for this purpose, are high- 
resolution Compton scattering^^'^- and two-dimensional 
angular correlation of positron annihilation radiation 
(2D-ACAR)ji2r— Compton scattering probes the mo- 
mentum density of the many-body electronic ground 
state of the system and is insensitive to the presence 
of defects or surfaces in the samplci^ 2D-ACAR also 
probes the bulk momentum density, but its interpreta- 
tion can be complicated by positron spatial distribution 
effects.""— On the other hand, the Compton scattering 
technique requires large single crystals of materials hav- 
ing atoms with high atomic numbers {Z ^ 30), as is the 
case for all cuprates, but it suffers from the problem of a 
relatively low signal from valence electrons sitting on the 
large background contribution from the core electrons. 

With this background, the present article reports a 



study of the FS and electron momentum density (EMD) 
of an overdoped single crystal of La2_2;Sr^Cu04 (LSCO) 
for the hole doping level x—0.3. Compton scattering 
and 2D- AGAR experiments have been carried out on the 
same LSCO sample, and the results are analyzed through 
parallel computations of the FS, the EMD as well as the 
electron-positron momentum density within the frame- 
work of DFT. The conventional picture of the metallic 
state based on Landau Fermi liquid theory is expected 
to become increasingly viable with doping as the system 
becomes more weakly correlated, even though the physics 
of the cuprates is generally dominated by deviations from 
such a simple picture of electronic states. For this reason, 
the overdoped LSCO provides a good starting system for 
investigating the fermiology of the cuprates. 

Although earlier Compton studies have provided some 
insight into the ground-state momentum density and 
electron correlation effects in LSCO^iiS^ we attempt to 
study the FS of LSCO using high-resolution Compton 
scattering experiments. The 2D-ACAR experiments, on 
the other hand, have been deployed successfully by sev- 
eral groups in the past for delineating the FSs of the 
cuprates in a number of favorable cases<^ 

Here, by measuring a series of high resolution Comp- 
ton profiles, we reconstruct the 2D momentum density in 
overdoped LSCO and identify a clear signature of the FS 
in the third Brillouin zone. Moreover, the DFT-based 
theoretical EMD is found to be in quantitative accord 
with the Compton scattering results, indicating that the 
ground state wavefunction in the overdoped system is 
well-described by the weakly correlated DFT picture. We 
have also found a quantitative level of agreement be- 
tween the measured and computed 2D-ACAR spectra. 
However, the 2D-ACAR spectra do not reveal clear FS 
signatures due to the well-known positron spatial distri- 
bution effects, which can make the positron insensitive 
to electrons in the Cu-0 planes. 

The remainder of this article is organized as fol- 
lows. Section H describes experimental details of sam- 
ple preparation and of Compton scattering and positron- 
annihilation experiments. Section HI provides technical 
details of our DFT-based EMD, FS and electron-positron 
momentum density computations. Section IV discusses 
momentum density anisotropy results, while Sec. V con- 
siders the FS determination. The article concludes with 
a summary of the results in Sec. VI. 



II. EXPERIMENTS 



evaporation of CuO during the high temperature pro- 
cess. The grown crystal was subsequently annealed un- 
der an oxygen pressure of 3 atm at 1173 K for 100 hours. 
Superconducting quantum interference device (SQUID) 
measurements, using MPMS-XL5HG (Quantum Design, 
Inc.), showed no superconductivity down to 2 K. Neutron 
diffraction studies indicated that the crystal is tetragonal 
(I4/mmm) down to the lowest temperature. The crystal 
was also characterized by other experiments i^"— 

We have measured 10 Compton profiles with scatter- 
ing vectors equally spaced between the [100] and [110] 
directions using the Cauchois-type x-ray spectrometer at 
the BL08W beamline of SPring-S.^'i All measurements 
were carried out at room temperature. The overall 
momentum resolution is estimated to be 0.13 a.u. 
full- width- at- half- maximum (FWHM). The incident 
x-ray energy was 115 keV and the scattering angle was 
2.88 rad. Approximately 5x10^ counts in total were 
collected at the Compton peak channel, and two inde- 
pendent measurements were performed in order to check 
the results. Each Compton profile was corrected for 
absorption, analyzer and detector efficiencies, scattering 
cross section, possible double scattering contributions, 
and x-ray background. The core-electron contribu- 
tions were subtracted from each Compton profile. A 
two-dimensional momentum density, representing a 
projection of the three-dimensional momentum density 
onto the a — b plane, was reconstructed from each set of 
ten Compton profiles using the direct Fourier transform 
methodiiSli 

The 2D-ACAR was measured using the Delft Univer- 
sity 2D-ACAR spectrometer^A with a conventional ^^Na 
positron source. The data were taken at a temperature 
of about T = 70 K. To correct for the sample shape 
(5.7 X 3.5 X 4.5 mm^) the data were convoluted with a 
gaussian of FWHM of 0.1 channel in the x direction and 
1.7 channels in the y direction (along which the positrons 
impinge on the sample), where one channel corresponds 
to 0.184 mrad. The total resolution is 1.0 x 1.0 mrad^ 
FWHM (1 mrad = 0.137 a.u). The total number of coin- 
cidences collected is 76.3 x 10^ and the maximum number 
of coincidences at the center is 16.7 x 10^ counts. 

The effects of superconductivity on the 2D-ACAR as 
well as Compton spectrai22 are expected to be below the 
resolution of the current experiments. For this reason, we 
have not carried out experiments on the superconducting 
state in connection with the present study. 



The heavily overdoped single crystal (x — 0.30) was 
grown by the traveling solvent fioating zone method. For 
this purpose, a powder sample was first synthesized by 
the conventional solid state reaction method. It was then 
shaped into feed rods under hydrostatic pressure and 
sintered at 1173 K for 12 hours, and at 1423 K for an 
additional 10 hours. In this process, excess CuO of 2 
mol% was added to the feed rods to compensate for the 



III. CALCULATIONS 

Our electronic structure calculations are based on the 
LDA within the framework of the DFT. An all-electron 
fully charge-self-consistent semi-relativistic Korringa- 
Kohn-Rostoker (KKR) methodology was used.— The 
crystal structure of LSCO was taken to be body-centered 
tetragonal (BCT) with space group lA/mmm (139) us- 



ing lattice parameters a and c given in Ref. |3J. A non- 
spinpolarized calculation neglecting the magnetic struc- 
ture was performed. Self-consistency was obtained for 
X — and the effect of doping was treated within a rigid 
band model by shifting the Fermi energy to accommo- 
date the proper number of electronsi^"— The results are 
in good agreement with other calculations.— The formal- 
ism for computing momentum density, p{p), is discussed 
in Refs. 139^421. The EMD is calculated according to the 
formula 



P(P) = E»"*(k) I J cxp{~ip.r)ij,{r)dr 



and the AGAR spectrum is computed as 

P^''{P) == Ej^i^k) I /exp(-ip.r)-0i(r)0+(r)dr |2, 



(1) 



(2) 



where ipi denotes the electronic wave function, 0+ the 
positron wavefunction, and p = k + G with G a recip- 
rocal lattice vector. ni(k) is the occupation function^ 
which in the independent particle model equals 1 if the 
electron state i is occupied and when it is empty, and 
the summation extends over the occupied k states. In 
the /O^^(p) calculation, we have neglected the enhance- 
ment factor for the annihilation rate.~ The inclusion of 
enhancement effects is crucial for the calculation of life- 
times but is well-known^^'*'^ to be not very important for 
discussing questions of bonding and FS signals in mo- 
mentum density. The momentum densities are calcu- 
lated on a momentum mesh with step {Spx,Spy,Spz) = 
27r(l/32a, l/32a, l/4c). The momentum is expressed in 
atomic units defined by 1/ao where ao is the Bohr ra- 
dius. The calculations include contributions from both 
the filled valence bands and the conduction-band which 
gives rise to the FS in LSCO. To study the electronic 
structure of the system, we consider two quantities of in- 
terest, 2D-EMD and ID-EMD, which are the projection 
of EMD in 2D and ID, respectively, given by: 



P^'^{Px,Py) = / p{p)dpz 



and 



rip.) 



j j pip)dpydpz 



(3) 



(4) 



The calculated band structure of LSCO {x = 0.3) near 
the Fermi level is illustrated in Fig. [IJa) . The band clos- 
est to the Fermi level is shown by the red dotted curve. 
This Cu02 energy band is dominated by copper-oxygen 
d^2_y2 —px,y orbitals. The present calculation for x=0.3 
predicts that the vHS is above the Fermi energy, so the 
FS has a diamond shape closed around the F point. Fig- 
ure mb) shows the 2D momentum density contribution 
of the x"^ —y"^ band, Px2 _y2 {Px,Py), together with the FS 
sections at (fcz = 0) for the doping level x = 0.3 mapped 
periodically throughout the momentum space. The Bril- 
louin zones are visualized in the figure by a grid of 27r/a, 
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FIG. 1: (Color online) (a) Band structure of LSCO near the 
Fermi level. The Cu02 band is shown as a red dotted line. 
(b) Calculated p^2_y2{Px,Py) is shown together with the FS 
sections at {kz — 0) for the doping level x — 0.3. The color 
scale is in units of p^''(0,0). The grid represents Brillouin 
zones with a size of 2n/a, where a is the lattice constant. 



where a is the lattice constant. The momentum density 
acts as a matrix element on the FS map. Thus, since 
liPxiPy) has strong amplitudes in the third Bril- 



Pit 



louin zones, the FS could be more easily detected there. 

The positron annihilation and Compton scattering ex- 
periments probe all the electrons in the system. However, 
core and semi-core electrons give isotropic distributions 
while the anisotropy of the spectra is produced by elec- 
trons near the Fermi energy. Therefore, we can concen- 
trate on analyzing the residual anisotropy after subtrac- 
tion of the isotropic part. We consider an anisotropy with 
C4y symmetry, given by^: 
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FIG. 2: (color online) Theoretical and experimental profiles FIG. 3: (color online) Theoretical and experimental y4^'*(p^ 
P^'^iPx) for (a) AGAR and (b) Compton scattering along the for (a) AGAR and (b) Gompton scattering 
[100] direction. All spectra are normalized to p^'^{0). 



In order to enhance FS signals, it is also possible to cal- 
culate another anisotropy obtained by subtraction of a 
smooth cylindrical average of the distribution p^'* {p^ , Py ) 
defined by 



A^\p.,Py) - p'^P.^Py) Si^pl+pl), 



(6) 



where 5* is a smoothed cylindrical average of p, in which 
the original spectrum is averaged over rotation angles 
from 0° to 45° in steps of 1°. 



IV. MOMENTUM DENSITY ANISOTROPY 

We start by comparing the experimental ID profiles 
based on Compton scattering and AGAR measurements 
with the corresponding theoretical predictions. Figure [2] 
shows a good level of accord between theory and exper- 
iment, especially for the Compton results. Further in- 
sight is gained by considering the ID-projection of the 
anisotropy of Eq. [S] given by 



A'^{P.) 



Pb 



AcuSP^^Pv)dPv 



(7) 



where Ap = Pb — Pa is the momentum range over which 
the projection is taken.4'^ The profile A^'^{px), shown in 
Fig. [21 is also equal to the difference of profiles^S between 
two crystallographic directions [100] - [110] and it can be 
compared to a similar Compton profile anisotropy mea- 
sured by Shukla et al. at a lower hole doping4^ The 
amplitude of the theory is the same as that in the ex- 
periment, while in Ref. |4^ the theoretical A^'''{px) had 
to be scaled down by a factor of 1.4 to obtain agreement 
with experiment. We note that theory predicts that the 
main contribution of correlation effects is an isotropic 
redistribution of the momentum density,"^- so that the 
amplitudes of oscillation in A^'^{px) become significantly 
reduced in the strongly correlated system. Hence our 
quantitative agreement between theory and experiment 



suggests that the x = 0.30 doped regime is consistent 
with Fermi liquid physics and that correlation effects 
modifying the anisotropy^^^ are relatively weak. 

In order to focus on the Cu-0 band contribution, 
Fig. m presents the 2D dy anisotropy distributions in 
the {px,Py) plane for positron annihilation and Comp- 
ton scattering spectra, respectively.^^ The anisotropy of 
the ACAR^ in Fig. 2] can be modeled by a molecular 
orbital method^^ involving the overlap of the positron 
wavefunction with Cu 3rf states hybridized with O 2p 
states. For an atomic orbital, the momentum density has 
the same point symmetry as the corresponding charge 
density. This result carries over to molecular states^ 
and is equally applicable to solid-state wave functions;^ 
The Compton scattering anisotropy maps are very sim- 
ilar overall to the positron-annihilation results, except 
that the Compton spectra extend to significantly higher 
momenta. This is expected since in the ACAR case the 
tendency of the positron to avoid positively charged ionic 
cores has the effect of suppressing higher momentum 
components of the EMD produced by the core and the 
localized valence electrons. Figure H] clearly shows that 
the theory reproduces most of the structure observed in 
the measured Compton as well as positron-annihilation 
distributions, including the momenta at which various 
features are located. The fourfold symmetry of the spec- 
tra is a consequence of the symmetry of the underlying 
body-centered tetragonal lattice. 

Therefore, an important finding of our work is that 
the correlations that ordinarily reduce the Compton 
anisotropy amplitudes are no longer effective at this dop- 
ing. Similar Fermi liquid behavior has been reported 
in studies of overdoped Tl-cuprates.J^^ Moreover, the 
trend of the weakening of correlation effects with dop- 
ing is also consistent with the changes observed in x-ray 
absorption^^'^^ and photoelectron spectroscopiesi^ 
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FIG. 4: (color online) Top: (left) Experimental and (right) 
theoretical AGAR C4v anisotropy. Bottom: (left) Experimen- 
tal and (right) theoretical C4v anisotropy (Compton). The 
color scale is in units of p^'*(0,0). 



V. FERMI SURFACE RESULTS 

We comment briefly on the positron results first. Al- 
though the electron-positron momentum density mea- 
sured in a positron-annihilation experiment contains FS 
signatures, the amplitude of such signatures is controlled 
by the extent to which the positron wavefunction over- 
laps with the states at the Fermi energy. Figure [S] shows 
the positron density distribution in a planar section of 
the LSCO unit cell. As in the calculation by Blandin 
et ali^'^, our result indicates that the positron does not 
probe well the FS contribution of the Cu-0 planes in 
LSCO. Moreover, positron-annihilation favors FS contri- 
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FIG. 5: (color online) Positron density for LSGO in the (010) 
plane. The positron wave-function is normalized in the LSGO 
unit cell and a logarithmic density scale (with integer numbers 
as units) is used in order to enhance regions with low positron 
density. The atomic positions (La, Gu, O) set the scale for 
the X and y axes. 



butions involving O 2p rather than Cu 3d states because 
the positron wavefunction overlaps more with extended 
O 2p states compared to the more tightly bound Cu 3d 
states. 

Indeed, we see little evidence of FS signatures in ei- 
ther the computed or the measured AGAR distributions 
FigO (a) . A more favorable case for the 2D-ACAR distri- 
bution is provided by the YBa2Cu307 cuprate supercon- 
ductor, where the 1-dimensional ridge FS has a two-fold 
symmetry which distinguishes it from important four- 
fold symmetry wave function effects j^^'^^ On the other 
hand, the Compton scattering calculation reveals a clear 
Fermi surface feature near a momentum (1.5, 1.5) a.u., 
and equivalent positions, Fig. 4(b). Strikingly, a similar 
feature is seen in the experimental spectrum on the left- 
side of the figure. To see this more clearly, we apply a 
Lock-Crisp- West folding procedure. Fig. 6. 

The Lock-Crisp- West (LCW) theorem— can be used 
both in Compton scattering 2D-EMD and in the case 
of 2D-ACAR data to study the non isotropic features of 
the momentum density by folding the data into a single 
central Brillouin zone. This technique can enhance FS 
discontinuities ('breaks') by coherently superposing the 
umklapp termsi^ However, the LCW folding can also ar- 
tificially enhance errors in the experimental dataJ^ Thus, 
to more clearly expose FS features in the data, we con- 
sider the cylindrical anisotropy defined by Eq. [S] Since 
the subtracted function is smooth and slowly varying this 
procedure does not contribute to, nor create, new struc- 
tures in the spectrum remaining after subtraction. 

Figure [S] (a) shows the cylindrical anisotropy of the 
theoretical 2D-EMD spectrum. The momentum density 
is represented in the extended zone scheme. Because the 
FS is periodic, a complete FS must exist in each Bril- 
louin zone, but with its intensity modulated by matrix 
element effects, as in Fig. 1(b). For a predominantly 
d-wave FS, the matrix element effects will strongly sup- 
press spectral weight near F, so the FS breaks are most 
clearly seen in higher Brillouin zones. These FS breaks 
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FIG. 6: (color online) (a) Cylindrical anisotropy for theoret- 
ical Compton scattering momentum density. The white lines 
define the Brillouin zones, while blue squares indicate a fam- 
ily of zones where the FS features are particularly strong, (b) 
The blue squared regions are isolated from the rest of the 
spectra and folded back to the central region. The color scale 
is in units of p^"* (0,0). 



matrix element effects are substantially circumvented. 
The resulting FS map is shown in the center of Figl6](b), 
and again on a larger scale in FiglT] (a). Applying the 
same procedure to other Brillouin zone regions produces 
similar results. For instance, the four regions along the 
diagonal neighboring the central region can also yield full 
FS information, but here it is superimposed on strong 
momentum density features. 

In Figure [7] we compare the theoretical FS obtained 
via the aforementioned 'partial folding' procedure to the 
corresponding experimental result. Figure [7] (a) shows 
that the partial backfolding technique produces the FS 
of correct size and topology. The same procedure when 
applied to the experimental spectrum yields results of 
Fig|7] (d). To compare with theory, we have convoluted 
the spectrum in FiglTja) with 0.07 a.u and 0.13 a.u. res- 
olutions, leading to the distributions of Figsl7](b) and (c) 
respectively. The latter is close to the actual experimen- 
tal resolution used in this work. 

By comparing the four frames in Fig. [71 one can see 
that the experimental spectrum in Fig. [7] (d) shows FS 
structures consistent with the LDA theory shown in 
Fig. [7] (a). Clearly the FS discontinuity along the nodal 
direction is smeared when including resolution broaden- 
ing as shown in FiglT] (b) and Fi^ (c ). The best over- 
all agreement between theory and experiment lies be- 
tween Fig. 13 (b) and Fi^ (c ). This is also confirmed 
by the cuts taken along the nodal direction plotted in 
Fig 13 The remaining discrepancy between experiment 
and theory could result from intrinsic or extrinsic inho- 
mogeneity effects such as the appearance of local ferro- 
magnetic clusters about concentrated regions of dopant 
atoms, ^ which have been neglected in the present simu- 
lations. Our conclusion is that the closed FS for x = 0.30 
predicted by our LDA calculation is consistent both with 
the present Compton scattering data and with surface- 
sensitive ARPES resuhs.^ 



VI. CONCLUSIONS 



appear superimposed on the momentum density in the 
form of discontinuities which can occur in any Brillouin 
zone. In FiglHl (a), the Fermi surface breaks are the re- 
gions where the contours run closely together so that the 
electron momentum density varies rapidly at these loca- 
tions. Figure IS] (a) shows the calculated Fermi breaks 
in several Brillouin zones. In particular, in the third 
zones framed with blue squares, the arc-like features are 
theoretically predicted FSs associated with Cu-0 planes. 
Due to the tetragonal symmetry, a rotation of the spec- 
trum by 7r/2 will generate symmetry-related regions (blue 
squares) with equivalent strong FS features. These re- 
gions are isolated in Figl6](b). By performing a 'partial 
folding', that is, folding only these regions back into the 
first Brillouin zone, we produce a full FS, where strong 



We have performed momentum density measurements 
on a high quality overdoped LSCO sample using both 
Compton scattering and 2D-ACAR. First principles 
calculations were also performed for the corresponding 
spectra. The quantitative agreement between the calcu- 
lations and the experiment for ACAR as well as EMD 
anisotropics suggests that x = 0.3 overdoped LSCO 
can be explained within the conventional Fermi-liquid 
theory. Nevertheless, a FS signal was only clearly 
observed by Compton scattering in the third Brillouin 
zone along [100]. Our FS analysis confirms previous 
ARPES FS measurements^ showing an electron-like FS 
in the overdoped regime. This validation is important 
since we provide via deep inelastic x-ray scattering 
experiments a truly bulk-sensitive image of momentum 
density maps of electrons near the Femi level. In general, 
this momentum density information is difficult to extract 




FIG. 7: Result of the 'partial folding' procedure: (a) Non 
convoluted Theory; (b) Theory convoluted with experimental 
resolution 0.07 a.u.. (c) Theory convoluted with experimental 
resolution 0.13 a.u.; (d) Experiment. The white color corre- 
sponds to a weight of unity (occupied) while the black color 
corresponds to zero (unoccupied). 



from ARPES experiments due to difficulties associated 
with matrix clement effects and the well-known surface 
sensitivity of ARPES. 
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FIG. 8: (color online) Cuts of the distributions in Fig.[7]along 
the nodal direction. The amplitudes of theory and experiment 
are compared on the same scale. Curves are offset with re- 
spect to one another for clarity. From the bottom to the top, 
the first curve is the non-convoluted theory, the second is the 
convoluted theory with a resolution 0.07 a.u., the third one 
is the convoluted theory with a resolution 0.13 a.u. and the 
fourth curve is the experimental data. 



P.E. Bisson, P. Descouts, A. Dupanloup, A. A. Manuel, E. 
Perreard, M. Peter and R. Sachot, Helv. Phys. Acta 55, 
110 (1982); J. Mader, S. Berko, H. Krakauer and A. Bansil, 
Phys. Rev. Letters 37, 1232(1976); J. C. Campuzano, L. 
C. Smedskjaer, R. Benedek, G. Jennings, and A. Bansil, 
Phys. Rev. B43, 2788(1991). 

^® L. P. Chan, K. G. Lynn and D. R. Harshman, Mod. Phys. 
Lett. B 6, 617 (1992). 

^^ R.N. West, in Positron Spectroscopy of Solids, A. Du- 
pasquier and A.P.Mills Jr. eds, (lOS Press, Amsterdam 
1995) pp. 75-143 

^° B. Barbiellini, A. Koizumi, P. E. Mijnarends, W. Al- Sawai, 
Hsin Lin, T. Nagao, K. Hirota, M. Itou, Y. Sakurai and A. 
Bansil, Phys.Rev. Lett. 102, 206402 (2009); Y. Sakurai, 
Y. Tanaka, A. Bansil, S. Kaprzyk, A. T. Stewart, Y. Na- 
gashima, T. Hyodo, S. Nanao, H. Kawata, and N. Shiotani, 
Phys. Rev. Letters 74, 2252(1995). 

^^ L. C. Smedskjaer, A. Bansil, U. Welp, Y. Fang and K. 
G. Bailey, Physica C 192, 259 (1992); A.A. Manuel, A. 
Shukla, L. Hoffmann, T. Jarlborg, B. Barbiellini, S. Mas- 
sidda, W. Sadowski, E. Walker, A. Erb and M. Peter, J. 
Phys. and Chem. of Solids 56, 1951 (1995); A.A. Manuel, 
Handbook on the Physics and Chemistry of Rare Earths 
Vol. 30, edited by K.A. Gschneider, Jr., L. Eyring and M.B. 
Maple, (Elsevier Science B.V. 2000); P. E. Mijnarends, A. 
C. Kruseman, A. van Veen, H. Schut and A. Bansil, J. 
Physics: Condens. Matter 10, 10383(1998). 

^^ P. Blandin, S. Massidda, B. Barbiellini, T. Jarlborg, P. 
Lerch, A. A. Manuel, L. Hoffmann, M. Gauthier, W. Sad- 
owski, E. Walker, M. Peter, Jaejun Yu and A. J. Freeman, 
Phys. Rev. B 46, 390 (1992). 

2=* R. H. Howell, P. A. Sterne, M. J. Fluss, J. H. Kaiser, K. 
Kitazawa and H. Kojima, Phys. Rev. B 49, 13127 (1994). 



J. Laukkanen, K. Hamalainen, S. Manninen, A. Shukla, T. 

Takahashi, K. Yamada, B. Barbiellini, S. Kaprzyk and A. 

Bansil, J. Phys. Chem. Sohds 62, 2249 (2001). 

Y. Sakurai, M. Itou, B. Barbiellini, P. E. Mijnarends, R. 

S. Markiewicz, S. Kaprzyk, J.-M. Gillet, S. Wakimoto, M. 

Fujita, S. Basak, Yung Jul Wang, W. Al-Sawai, H. Lin, A. 

Bansil and K. Yamada, Science 332, 698 (2011). 

S. Wakimoto, H. Zhang, K. Yamada, I. Swainson, 

Hyunkyung Kim and R. J. Birgeneau, Phys. Rev. Lett. 

92, 217004 (2004). 

S. Wakimoto, Young- June Kim, Hyunkyung Kim, H. 

Zhang, T. Gog and R. J. Birgeneau, Phys. Rev. B 72, 

224508 (2005). 

S. Wakimoto, R. J. Birgeneau, A. Kagedan, Hyunkyung 

Kim, I. Swainson, K. Yamada and H. Zhang, Phys. Rev. 

B 72, 064521 (2005). 

N. Hiraoka, M. Itou, T. Ohata, M. Mizumaki, Y. Sakurai 

and N. Sakai, J. Synchrotron Radiat. 8, 26 (2001); M. 

Itou and Y. Sakurai, AIP Conf. Proc. 705, 901 (2004); 

Y. Sakurai and M. Itou, J. Phys. Chem. Solids 65, 2061 

(2004). 

1. Matsumoto, J. Kwiatkowska, F. Maniawski, M. Itou, H 
Kawata, N. Shiotani, S. Kaprzyk, P. E. Mijnarends, B 
Barbiellini and A. Bansil, Phys. Rev. B 64, 045121 (2001) 
C. V. Falub, P. E. Mijnarends, S.W.H. Eijt, M.A. van Huis 
A. van Veen and H. Schut, Phys. Rev. B 66, 075426 (2002) 
S.E. Barnes and M. Peter, Phys. Rev. B 40, 10958 (1989) 
A. Bansil, S. Kaprzyk, P. E. Mijnarends and J. Tobola 
Phys. Rev. B 60, 13396 (1999). 

S. Sahrakorpi, M. Lindroos, R. S. Markiewicz and A. Ban- 
sil, Phys. Rev. Lett. 95, 157601 (2005). 
The rigid band picture^ is expected to be reasonable for 
doping away from the cuprate planes. It will be interesting 
to examine doping effects via first principles approaches^. 
A. Bansil, Z. Naturforsch. 48 A, 165 (1993); Hsin Lin, S. 
Sahrakorpi, R. S. Markiewicz and A. Bansil, Phys. Rev. 
Letters 96, 097001 (2006); H. Asonen, M. Lindroos, M. 
Pessa, R. Prasad, R. S. Rao and A. Bansil, Phys. Rev. 
B25, 7075 (1982). 

S. Kaprzyk and A. Bansil, Phys. Rev. B 42, 7358(1990); 
L. Schwartz and A. Bansil, Phys. Rev. B 10, 3261(1974); 
R. Prasad and A. Bansil, Phys. Rev. B21, 496 (1980). 
J. Yu, A.J. Freeman and J.-H. Xu, Phys. Rev. Lett. 58, 
1035 (2008). 

P. E. Mijnarends and A. Bansil, Phys. Rev. B 13, 2381 
(1976); ibid. Phys. Rev. B 19, 2912 (1979). 
P. E. Mijnarends and A. Bansil, J. Phys.: Condens. Matter 

2, 911 (1990). 

A. Bansil, R. S. Rao, P. E. Mijnarends and L. Schwartz, 
Phys. Rev. B 23, 3608 (1981). 

P.E. Mijnarends and A. Bansil, in Positron Spectroscopy 
of Solids, edited by A. Dupasquier and A. P. Mills (lOS 
Press, Amsterdam, 1995), p. 25. 

B. Barbiellini and A. Bansil, J. Phys. Chem. Solids 62, 
2181 (2001); B. Barbiellini, J. Phys Chem. Solids 61, 341 
(2000). 

B. Barbiellini, M. J. Puska, T. Torsti and R. M. Nieminen, 
Phys. Rev. B 51, 7341 (1995); B. Barbiellini, M. J. Puska, 
T. Korhonen, A. Harju, T. Torsti and R. M. Nieminen, 
Phys. Rev. B. 53, 16201 (1996). 

B. Barbiellini, P. Genoud, P. Lerch, T. Jarlborg and M. 
Peter, Materials Sciences Forum 105-110, 599 (1992). 
C4v is the Schoenflies notation for a group containing a 
a„ refiection in addition to C4, a single 4- fold axis sym- 



Manuel, T. Bus- 
Rev. B 59, 12127 

B 9, 5122 (1974). 
Sci. Forum 255- 



metry. See e.g. M. Tinkham, Group Theory and Quantum 
Mechanics, (McGraw-Hill, 1964). 

The limits of integration are within the momentum interval 
[-2.5, 2.5] a.u. 

A. Shukla, B. Barbiellini, A. Erb, A.A. 
laps, V. Honkimaki and P. Suortti, Phys. 
(1999). 
L. Lam and P.M. Platzman, Phys. Rev. 

B. Barbiellini and P.M. Platzman, Mat 
257, 189 (1997). 

Our experimental AGAR spectra are in reasonable accord 
with earlier studies Refs. I22ll23l . 

P.E.A. Turchi, A.L. Wachs, K.L. Wetzler, J.H. Kaiser, 
R.N. West, Y.G. Jean, R.H. Howell and M.J. Fluss, J. 
Phys. Gondens. Matter 2, 1635 (1990). 
R. Harthoorn and P.E. Mijnarends, J. Phys. F 8, 1147 
(1978). 

T. Ahmed, T. Das, J.J. Kas, Hsin Lin, B. Barbiellini, Fer- 
nando D. Vila, R.S. Markiewicz, A. Bansil and J. J. Rehr, 
Phys. Rev. B 83 115117 (2011). 

D. G. Peets, D. G. Hawthorn, K. M. Shen, Y.-J. Kim, D. 
S. Ellis, H. Zhang, S. Komiya, Y. Ando, G. A. Sawatzky, 
R. Liang, D. A. Bonn and W. N. Hardy, Phys. Rev. Lett. 



103, 087402 (2009). 

^^ M. Schneider, R.-S. Unger, R. Mitdank, R. MiiUer, A. 
Krapf, S. Rogaschewski, H. Dwelk, G. Janowitz and R. 
Manzke, Phys. Rev. B 72, 014504 (2005). 

^^ H. Haghighi, J. H. Kaiser, S. Rayner, R. N. West, J. Z. 
Liu, R. Shelton, R. H. Howell, F. Solal, P. A. Sterne and 
M. J. Fluss, Phys. Rev. Lett. 67, 382 (1991). 

^^ L. Hoffmann, A.A. Manuel, M. Peter, E. Walker, M. Gau- 
thier, A. Shukla, B. Barbiellini, S. Massidda, Gh. Adam, S. 
Adam, W. N. Hardy and Ruixing Liang, Phys. Rev. Lett. 
71, 4047 (1993). 

^^ D.G. Lock, V.H.G. Grisp and R.N. West, J. Phys. F 3, 561 
(1973). 

^° T. Yoshida, X.J. Zhou, M. Nakamura, S.A. Kellar, P.V. 
Bogdanov, E.D. Lu, A. Lanzara, Z. Hussain, A. Ino, T. 
Mizokawa, A. Fujimori, H. Eisaki, G. Kim, Z.X. Shen, T. 
Kakeshita and S. Uchida, Phys. Rev. B. 63, 220501 (2001); 
T. Yoshida, X. J. Zhou, K. Tanaka, W. L. Yang, Z. Hus- 
sain, Z. X. Shen, A. Fujimori, S. Sahrakorpi, M. Lindroos, 
R. S. Markiewicz, A. Bansil, Seiki Komiya, Yoichi Ando, 
H. Eisaki, T. Kakeshita and S. Uchida, Phys. Rev. B 74, 
224510 (2006). 



